Planets in Stellar Clusters Extensive Search. II. Discovery of 57 
Variables in the Cluster NGC 2158 with Millimagnitude Image 

Subtraction Photometry. 
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ABSTRACT 

We have undertaken a long-term project, Planets in Stellar Clusters Extensive Search 
(PISCES), to search for transiting planets in open clusters. NGC 2158 is one of the targets 
we have chosen - an intermediate age, populous, rather metal poor cluster. In this paper we 
present the results of a search for variable stars in the data from the first season of monitoring at 
the FLWO 1.2 m telescope. This is the first variability search ever conducted in this cluster. We 
present a catalog of 57 variable stars, most with low amplitude variability. Among the variables is 
a cataclysmic variable (CV) which underwent a 2.5 mag outburst. If it is a member of NGC 2158, 
this would be the fourth CV known in an open cluster. We have also found five 5 Scuti stars, 
three of which we have two or more detectable modes of pulsation. Of the 57 variables discovered, 
28 have i?-band amplitudes of 5% or below. Six of those vary at or below the 2% level, including 
one with 0.08% variability. 
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We have undertaken a long-term project, Plan- 
ets in Stellar Clusters Extensive Search (PISCES), 
to search for transiting planets in open clusters. 
In Mochejska et al. (2002, hereafter Paper I) we 
demonstrated that we are capable of attaining the 
requisite photometric precision, based on one sea- 
son of data for NGC 6791. 

In this paper we present the results of a search 
for variable stars in our second target, NGC 2158, 
based on the data from the first observing season. 

NGC 2158 is a rich open cluster of intermediate 
age (2-3 Gyr). It is rather metal poor ([Fe/H]=- 
0.46) and is located at a distance of 3600 pc (Car- 
raro et al. 2002, hereafter Ca02; Christian, Heasley 
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and Janes 1985). Despite being one of the most 
populous galactic open clusters, it has never been 
investigated for variability. The only published 
CCD photometry for this cluster is that of Chris- 
tian et al. (1985), Piersimoni et al. (1993) and 
Ca02. An extensive catalog of photographic pho- 
tometry and proper motions has been published 
by Kharchenko, Andruk & Schilbach (1997, here- 
after KAS97). 

The paper is arranged as follows: §2 describes 
the observations, §3 summarizes the data reduc- 
tion procedure and §4 contains the variable star 
catalog. Concluding remarks are found in §5. 

2. Observations 

The data were obtained at the Fred Lawrence 
Whipple Observatory (FLWO) 1.2 m telescope us- 
ing the 4Shooter CCD mosaic with four thinned, 
back side illuminated AR coated Loral 2048 2 
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Fig. 1 — Digital Sky Survey image of NGC 2158 
showing the field of view of the 4Shooter. The 
chips are numbered clockwise from 1 to 4 starting 
from the bottom left chip. NGC 2158 is centered 
on Chip 3. North is up and east is to the left. The 
fields covered by KAS97 and Ca02 are plotted with 
dashed and dot-dashed lines, respectively. 



CCDs (Szentgyorgyi et al. 2004). The camera, 
with a pixel scale of 0'.'33 pixel -1 , gives a field 
of view of 1 1 f 4 x 11. '4 for each chip. The clus- 
ter was centered on Chip 3 (Fig. 1). The data 
were collected during 20 nights, from 3 January 
to 11 March, 2003. A total of 250 x 900 s R and 
69 x 450 s T^-band exposures were obtained. The 
median seeing was 1'.'9 in R and 2'.'1 in V. 

3. Data Reduction 

3.1. Image Subtraction Photometry 

The preliminary processing of the CCD frames 
was performed with the standard routines in the 
IRAF ccdproc package. 1 

Photometry was extracted using the ISIS image 
subtraction package (Alard & Lupton 1998, Alard 
2000). A brief outline of the applied reduction 
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procedure is presented here. For a more detailed 
description the reader is referred to Paper I. 

The ISIS reduction procedure consists of the 
following steps: (1) transformation of all frames 
to a common (a;, y) coordinate grid; (2) construc- 
tion of a reference image from several of the best 
exposures; (3) subtraction of each frame from the 
reference image; (4) selection of stars to be pho- 
tometered and (5) extraction of profile photometry 
from the subtracted images. 

All computations were performed with the 
frames internally subdivided into four sections 
(sub_x=sub_y=2). Differential brightness varia- 
tions of the background were fit with a first de- 
gree polynomial (deg_bg=l). A convolution ker- 
nel varying quadratically with position was used 
(deg_spatial=2). The psf width (psf _width) 
was set to 33 pixels and the photometric radius 
(radphot) to 5 pixels. The reference images were 
constructed from 20 best exposures. 

3.2. Calibration 

The VR photometry in Chips 3 and 4 was cal- 
ibrated from the photometry published by Ca02. 
1410 stars above R = 18.5 and 1312 above V = 19 
were used to determine the zero point of the mag- 
nitude scale for Chip 3. Since the overlap with 
Chip 4 was marginal (see Fig. 1), we used 44 stars 
above R = 19 and V — 20. The rms scatter of the 
residuals in R and V, respectively, was 0.046 and 
0.052 in Chip 3 and 0.065 and 0.059 in Chip 4. 

The i?-band photometry for Chips 1 and 2 was 
calibrated from the KAS97 photographic photom- 
etry. To determine the zero point we used stars 
above 17 mag in both bands: 39 and 19 in R and 
V for Chip 1, and 169 and 87 for Chip 2. The rms 
scatter of the residuals in R and V, respectively, 
was 0.18 and 0.14 in Chip 1 and 0.17 and 0.12 in 
Chip 2. 

As a consistency check we have calibrated Chip 
3 and 4 photometry from the KAS97 catalog and 
found an offset between the calibration based on 
Ca02. We have also compared the KAS97 and 
Ca02 catalogs directly, and confirm the difference 
between their zero points. The offsets range from 
0.11 to 0.19 in R and 0.04 and 0.1 in V, depend- 
ing on the sample of stars used. Since we could 
not constrain well the offsets between KAS97 and 
Ca02, we have left Chips 3 and 4 on the Ce02 scale, 



2 




Fig. 2. — V/V — R color-magnitude diagram (CMD) for Chip 3, centered on NGC 2158. Eclipsing binaries 
are plotted with circles, other periodic variables with squares and miscellaneous variables with triangles. 
The CV candidate, V57, is plotted at maximum and minimum (filled and open symbol, respectively). There 
are 3842 stars above R = 21. 



and Chips 1 and 2 on the KAS97 scale. 

In Figures 2 and 3 we present the V/V — R 
color-magnitude diagrams (CMD) for Chip 3 and 
Chips 1, 2, 4, respectively. The CMD for Chip 3, 
centered on the cluster, shows a well defined main 
sequence (MS) from the turnoff down to V ~ 21. 
The CMDs for the remaining chips are dominated 
by disk stars. There are 1429, 1764, 3842 and 1436 
stars above R = 21 in Chips 1 through 4, respec- 
tively. 

3.3. Astrometry 

Equatorial coordinates were determined for the 
i?-band template star lists. The transformation 
from rectangular to equatorial coordinates was de- 
rived using 740, 922, 1170, and 777 transformation 
stars from the USNO B-l catalog (Monet et al. 



2003) in Chips 1 through 4, respectively. The av- 
erage difference between the catalog and the com- 
puted coordinates for the transformation stars was 
0'/09 - O'.'ll. 

3.4. Selection of variables 

We extracted the light curves for all stars de- 
tected by DAOphot (Stetson 1987) on the tem- 
plate frames and searched them for variability us- 
ing the index J (Stetson 1996), as described in 
more detail in Paper I. The definition of J requires 
observations to be closely spaced in time. Points 
separated by a value smaller than some threshold 
(0.03 days in our case) are grouped into pairs and 
only when both have a residual from the mean of 
the same sign, the pair contributes positively to 
the variability index. 
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Fig. 3.— V/V — R CMDs for Chips 1, 2 and 4. Eclipsing binaries are plotted with circles, other periodic 
variables with squares and miscellaneous variables with triangles. There arc 1429, 1764 and 1436 stars above 
R = 21 in Chips 1, 2 and 4, respectively. 



To search for periodicities we have used the 
method introduced by Schwarzenberg-Czcrny 
(1996), employing periodic orthogonal polynomi- 
als to fit the observations and the analysis of vari- 
ance (ANOVA) statistic to evaluate the quality of 
the fit. 

4. Variable Star Catalog 

We have discovered 57 variables in all four 
chips, all newly discovered: 34 eclipsing binaries, 
19 other periodic and 4 long period or non-periodic 
variables. 

The properties of the variables are discussed 
in §4.1-4.5. Tables 1, 2 and 3 list the follow- 
ing parameters for the variables: the identifica- 
tion number, right ascension, declination, period 
(if applicable), R and V-band magnitudes (maxi- 
mum for non-periodic and eclipsing variables, flux 
weighted mean for periodic) and variability ampli- 
tudes in each band (semiamplitudes for periodic 
variables). The last two columns in Tab. 1 and 
2 provide the probabilities Pi and Pi for cluster 
core and corona membership, taken from KAS97. 
The i?-band light curves of all variables are shown 
in Figures 4- 8. 2 

The variables are plotted on CMDs for their re- 

2 The VR band photometry and finding charts for all vari- 
ables are available from the authors via the anonymous ftp 
on cfa-ftp.harvard.edu, in the /pub/bmochejs/PISCES di- 
rectory. 



spective chips (Figs 2 and 3). Eclipsing binaries 
are marked with circles, other periodic variables 
with squares and miscellaneous variables with tri- 
angles. By far the largest number of variables, 
34, were found on Chip 3, centered on the cluster, 
compared to 6, 7 and 10 on Chips 1, 2 and 4, re- 
spectively. This is in agreement with the 5' radius 
estimate for the cluster (Lynga 1987). 

4.1. Eclipsing Binaries 

We have identified 34 eclipsing binaries. Six 
are contact W UMa type systems (V01, V02, V04, 
V05, V07, V12). More than half of the variables 
could either be EB-type eclipsing binaries or ellip- 
soidal variables, as most have periods in excess of 
a day and amplitudes below 0.1 mag (V03, V08, 
VI 1, V14-V27, V33). It is possible that some of 
the longer period systems could be BY Dra type 
spotted stars with half the reported period. There 
are also seven Algol type detached binaries (V09, 
V10, V13, V28-V30, V32, V34). V06 and proba- 
bly V31 are RS CVn type binaries, where one of 
the components has non-uniform surface bright- 
ness due to cool spots. A starspot wave is seen on 
top of the variability due to eclipses. 

Proper motion data from KAS97 show that 
V12, V22 and V24 are most likely cluster mem- 
bers, while V14 is rather not. For V21, V28 and 
V31 the data are inconclusive. 
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Fig. 4. — The light curves of the 34 eclipsing binaries. 
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Fig. 4. — Continued. 



V35 



15.79 


; 1 1 


i i i i i i 

P=o!0981 




; ■. p=o!o98i • •. : 
.' ; .-■ sub :. •■; ; ; .: 


15.8 










15.81 










15.82 








7 ' ". * . ■ * ~ 


15.83 










15.79 
15.8 




1 'l 1 | 1 1 1 1 | 1 
P=o!l044 




1 i i 'i i i i i i | i i'l 

■ .. P=0!l044 .. 

: SUB . ; •'' : 


15.81 
15.82 




•V: ••V : "* '••'i'"' v. 




'- *•'.:■;■/.■ '■■lf*'- :r -6«*- 


15.83 


: I i 


, , , i ■, . , , i 




-i . , , , i ,• , , , i , , ,- 







0.5 1 




0.5 1 



phase 

Fig. 5. — The i?-band light curves of S Scuti vari- 
able V35, phased with each of the two detected 
periods: raw in the left panels, and subtracted 
of variability corresponding to other period in the 
right panels. 



4.2. 5 Scuti stars 

We have identified five variables, V35-V39, with 
periods 0.0981-0.1220 d and amplitudes 0.008- 
0.076 mag, located in the blue straggler region 
of the CMD. Their periods and amplitudes are 
typical for <5 Scuti variables. Their positions on 
the CMD are also consistent with their identifi- 
cation as S Scuti variables. In clusters consider- 
ably older than 1 Gyr the main sequence turnoff 
is located redward of the instability strip. Con- 
sequently, only blue stragglers can display <5 Scuti 
variability (Rodriguez & Breger 2001). 

For three of these variables, V35-V37, we have 
detected pulsations in more than one mode. Ta- 
ble 4 lists identified periods. Figures 5 and 6 show 
their light curves, phased with each detected pe- 
riod: raw in the left panels, and subtracted of vari- 
ability corresponding to other periods in the right 
panels. The simultaneous excitation of multiple 
modes of pulsation is typical for S Scuti variables, 
and such stars have been detected in open clus- 
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Fig. 6. — The i?-band light curves of 5 Scuti variables V36 and V37, phased with each of the three detected 
periods: raw in the left panels, and subtracted of variability corresponding to other periods in the right 
panels. 



ters (Rodriguez 2002, Frcyhammer, Arentoft & 
Sterken 2001, Pigulski, Kolaczkowski, & Kopacki 
2000). 

The ratios of the periods range from 1.020 to 
1.071. Such close frequency pairs are observed in 
the majority of well studied 5 Scuti stars (Breger 
& Bischof 2002). For periods around 0.1 d the 
ratio between the fundamental and first overtone 
period for Population I 5 Scuti stars is generally 
below 0.775 (Templeton, Basu & Demarque 2002). 
This implies that at least one of the excited modes 
of pulsation is non-radial. 

Using the fundamental mode Period-Luminosity 
(P-L) relation for 8 Scuti stars (Eq. 4 in Petersen 
& Christensen-Dalsgaard 1999), Rv = 3.2 (Patri- 
archy Morbidclli, & Pcrinotto 2003), E(B-V) = 
0.60 3 and (m - M) = 12.80 (Ca02), we get devi- 
ations from the cluster distance modulus less than 
0.15 mag for all stars except V39 (Tab. 4). This 



3 The final determination of E(B — V) in Ca02 was 0.55 ± 
0.10 mag. Here we are using the reddening and distance 
modulus obtained by them from isochrone fits. 



confirms that these stars are most likely 5 Scutis 
belonging to NGC 2158. One possible cause of 
discrepancy for V39 could be that the mode of 
pulsations we have identified is not the fundamen- 
tal mode. 

Four of these stars have proper motion data in 
KAS97. Cluster membership of V36 is confirmed 
with high probability, while for V35 and V39 the 
data are not conclusive. The catalog reports a very 
low cluster membership probability for V38, which 
is in disagreement with our conclusions based on 
the 6 Scuti P-L relation. The study of KAS97 was 
also targeting much brighter stars (V > 8 mag) in 
M 35 and two other clusters at smaller distances 
than NGC 2158. The authors say that their proper 
motion measurement "accuracy decreases rapidly 
for stars fainter than V — 15.5 mag", while the V 
magnitude range of NGC 2158 variables is 15.7- 
18.0, and V = 16.04 for V38. 
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4.3. Other periodic variables 

All of the remaining periodic variables are lo- 
cated redward of the MS, except for V42 and V43, 
which are located at its blue edge. Some of these 
variables could be BY Draconis type variables - 
rotating spotted stars. Six of the 13 variables are 
located on Chips 1, 2 and 4, making it likely that 
many of them are not cluster members. 

V48 seems to have a variable amplitude. Vari- 
ables V44, V48 and V49 on Chip 3, if they be- 
long to the cluster, might be members of the newly 
proposed class of variable stars termed "red strag- 
glers" (Albrow et al. 2001) or "sub-subgiant stars" 
(Mathieu et al. 2003). KAS97 gives a probability 
of corona memberships of 80%, 69% and 41% for 
V44, V48 and V49, respectively. To date, six such 
stars have been found in 47 Tuc (Albrow et al. 
2001) and two in M67 (Mathieu et al. 2003). Thus 
far, the origin and evolutionary status of these 
stars remains unknown. 

Variables V50 and V52 on Chip 1 and the RS 
CVn eclipsing binary V31 on Chip 4 also fall in 
the same region in the CMD. Since they are fur- 
ther from the cluster center, it is less likely they 
are members of NGC 2158 and of this class. The 
KAS97 catalog suggests V50 is most likely not a 
cluster member, and it is inconclusive for V31. 

The variable V42 has a sharply peaked light 
curve, reminiscent of RR Lyrae or short-period 
Cepheids. Its amplitude of 0.14 in R is much too 
low for either of those two types of variables, un- 
less it is blended with another star along the same 
line of sight. The KAS97 catalog gives it a high 
probability of being a cluster member. If V42 were 
really a member, it would have My ~ 3.2 and 
would be too faint to be either of these variables. 
On the CMD it is located at the blue edge of the 
MS (Fig. 2). The data for V42 from all nights, 
except the sixth and seventh, phase well with a 
period of 1.5393 days. In Fig. 7 the points from 
these two discrepant nights are marked with solid 
gray squares. For the full dataset (or even nights 
1-7) a good period cannot be determined. The ex- 
amination of images from nights six and seven did 
not reveal any defects, such as bad columns, ex- 
tremely high background or a bad psf, or anything 
that would suggest that the brightness measure- 
ments might be less accurate. Also, the data from 
these two nights show a coherent rise in brightness, 



which makes it more unlikely that these measure- 
ments are erroneous. 

4.4. Cataclysmic Variable Candidate 

We have identified a cataclysmic variable (CV) 
candidate V57. The star has undergone an out- 
burst with an amplitude of 2.5 mag in R and V. 
Its very blue color, V — R = —0.26 in maximum 
and —0.31 in minimum, supports this classifica- 
tion. If confirmed spectroscopically as a CV, and 
MS <X member of NGC 2158, this would be only 
the fourth CV known in an open cluster, with the 
other three being star 186 in M67 (Gilliland et al. 
1991) and B7 and B8 in NGC 6791 (Kaluzny et 
al. 1997). 

The V — R color of this CV candidate is in fact 
implausibly blue. We searched for this star in 
the Ca02 catalog, to check its color, but it is not 
listed there. We have obtained independent color 
derivations for this variable from three sets of R 
and V-band images taken on different nights when 
this object was near maximum brightness. One of 
these sets came from our most recent observations, 
not included in this analysis. These sets were cal- 
ibrated independently against the Ca02 catalog. 
The V — R color estimates ranged from -0.24 to - 
0.31. Unfortunately the seeing in V was over 2" on 
all V-band images when the variable was at max- 
imum, so the V-band magnitude probably has a 
considerable error associated with it. 

Assuming the NGC 2158 distance and redden- 
ing derived by Ca02, we get the minimum ab- 
solute V-band magnitude Mv(min) of 6.4 mag. 
From Fig. 3.5 in Warner (1995) it is apparent that 
My (min) is roughly 8 mag for U Gem type dwarf 
novae, but even the most extreme Z Cam type sys- 
tems have My (min) below 6.8 mag. This would 
imply that V57 is a dwarf nova located in front of 
the cluster, and would also make the color less ex- 
treme, assuming it has a much smaller reddening. 

It should be noted, however, that the V band 
magnitude alone is not sufficient to exclude this 
variable as a cluster member. Another cataclysmic 
variable, B7 in NGC 6791, is seen most of the time 
in a high state at My ~ 5, and has been caught 
only once in its low state at My (min) = 7.4. Fur- 
ther investigation is needed to resolve the nature 
of V57. A spectrum of a dwarf nova in quies- 
cence should display Balmer emission on a blue 
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continuum. During outburst the emission lines are 
gradually overwhelmed by the increasing contin- 
uum and development of broad absorption lines 
(Warner 1995). If V57 is in a "high" level, like 
B7, then this will be apparent in the spectrum 
and emission lines should not be prominent. 

4.5. Long period or non-periodic variables 

The variable V54 exhibits irregular variability 
and is located in the blue straggler region of the 
CMD. V55 is a detached eclipsing binary with one 
primary and one secondary eclipse observed, also 
probably a blue straggler. 

The shape of the light curve of V56 is very sim- 
ilar to V57, only with a much smaller amplitude 
of 0.07 in R and 0.05 in V. On the CMD it is 
located near the main sequence turnoff. We have 
considered the possibility that the variability of 
V56 could be an artifact related to V57, although 
the fact that V56 was detected on a different chip 
than V57 makes it less likely. We have examined 
images from two nights, when the variable was 
at minimum and maximum brightness, and found 
nothing suspicious. A simple examination under 
IRAF with imexam confirms the variability of this 



star, as well as its amplitude. 

5. Conclusions 

In this paper we present VR photometry for 
57 newly discovered variable stars in the inter- 
mediate age, populous, rather metal poor cluster 
NGC 2158. We have identified 34 eclipsing bina- 
ries, 19 other periodic variables and 4 long period 
or non-periodic variables. Among these variables 
we have discovered five 5 Scuti stars. For three 
of them we have identified two or more modes of 
pulsation. We have also found a cataclysmic vari- 
able which underwent a 2.5 mag outburst. If it is 
a member of NGC 2158 (which at this point does 
not seem very likely) this would be the fourth CV 
known in an open cluster. 

Of the 57 variables discovered, 28 have i?-band 
amplitudes of 5% or below. Six of those vary at 
or below the 2% level, including one with 0.08% 
variability. This, together with the results from 
Paper I, shows that we are capable of attaining 
requisite photometric precision to detect variabil- 
ity at the level expected for planetary transits. 
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Table 1 

Eclipsing binaries in NGC 2158 
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11 


47 
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V34 


6 
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23 
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0.040 


0.033 







Note. — Cluster core and corona membership probabilities P\ and Pi taken from K97. 
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Table 2 

Other periodic variables in NGC 2158 
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Note. — Cluster core and corona membership probabilities Pi and P2 taken from K97. 



Table 3 

Miscellaneous variables in NGC 2158 



ID 


02000 [h] 


*2000 [°] 


Rmax 
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6 07 07.0 


24 05 25.4 


15.337 


15.859 


0.034 


0.033 


V55 


6 07 25.7 


24 05 45.7 


15.630 


16.072 


0.370 


0.340 


V56 


6 06 17.2 


24 03 39.8 


16.699 


17.278 


0.074 


0.052 


V57 


6 07 33.8 


24 07 55.2 


18.847 


18.590 
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Table 4 

Properties of S Scuti stars in NGC 2158 



ID 


Pi [d] 


P 2 [d] 


P 3 [d] 


(m-M)o 


V35 


0.098051 


0.104419 




12.65 


V36 


0.107830 


0.111341 


0.114875 


12.89 


V37 


0.111285 


0.106032 


0.113531 


12.71 


V38 


0.120462 






12.66 


V39 


0.121997 






12.31 
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